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In swelling porous media, the potential for flow is much more than pressure, 
and derivations for flow equations have yielded a variety of equations. In this 
paper we show that the macroscopic flow potentials are the electro-chemical 
p • potentials of the components of the fluid and that other forms of flow equations, 

such as those derived through mixture theory or homogenization, are a result 
of particular forms of the chemical potentials of the species. It is also shown 
that depending upon whether one is considering the pressure of a liquid in a 
reservoir in electro-chemical equilibrium with the swelling porous media, or the 
pressure of the vicinal liquid within the swelling porous media, a critical pressure 
gradient threshold exists or does not. 
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1 Introduction 



Swelling porous materials are ubiquitous - they occur in soils such as swelling 
clays (montmorillonite), biotissues (cartilage), and in drug delivery systems such 
as Aleve (swelling polymers). Experiments are performed at the microscale 
(scale at which the solid and liquid or adsorbed liquid can be distinguished) 
and at the macroscale (scale at which the swelling porous media appears to be 
homogeneous, i.e. one cannot distinguish between the phases). The concept of 
pressure at each of these scales are often confused and interchanged. Example 
of terms used include 'disjoining pressure', 'osmotic pressure', and 'swelling 
pressure' and are attributed to the double-layer forces, van der Waals dispersion 
forces, osmotic forces, and surface hydration forces. 
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In addition, it is unclear which microscopic forces are dominant to 
macroscopic behavior. For example, although direct measurements indicate 
surface hydration is considered to have short-range effects (up to 4 monolayers 
of water entering between layers of montmorillonite clay) [34, 12], experiments 
by Low [28, 27] indicate that the macroscopic affects of these interactions can 
explain the osmotic swelling of montmorillonite soils in which the swelling is due 
to lOO's of layers of water. 

Several upscaling approaches have been used to arrive at a macroscopic 
model for flow through a swelling porous media, and with these approaches 
a variety of definitions of macroscopic flow potentials. It is the purpose 
of this paper to propose a macroscopic form for flow, derived from a 
hybrid mixture theory formulation [8, 9], and demonstrate how the form 
involving electrochemical potentials is a generalization of equations derived 
using homogenization [3(j], and is consistent with a Lagrangian mixture theoretic 
approach [23, 20]. In the process we illustrate that a pressure gradient threshold 
may exist, depending upon how the pressure is measured. 

For simplicity we assume that the swelling porous medium is composed of 
a solid and liquid phase (i.e. no gaseous phase). The solid phase (polymer, 
montmorillonite) is assumed to be negatively charged and the fluid contains 
cations, ions, and a neutral liquid. We will refer to the liquid phase as vicinal 
fluid to distinguish it from the bulk phase (liquid unaffected by its vicinity to 
the solid phase, or reservoir fluid). 

In the first section we review the microscale forces. We next review 
macroscopic quantities: osmotic repulsion, surface hydration, and disjoining 
pressure. We derive the flow equation in terms of chemical potentials from 
hybrid mixture theory results and discuss pressure gradient thresholds. Then we 
illustrate how the potential form of the flow equation can be used to derive forms 
derived via homogenization (Moyne and Murad [30]) used to model swelling 
montmorillonite, and is consistent with the mixture theory approach of Huyghe 
and Janssen [23] used to model swelling biotissues [20]. Although the models 
appear quite different, we show they can be derived from the potential form 
under particular assumptions on the chemical potentials of the species. 

2 Microscale Forces 

At the microscale there are various forces, some attractive and some repulsive, 
that cause a swelling porous medium to swell (repulsive forces dominate) or 
shrink (attractive forces dominate). In this section we summarize some of the 
forces considered to be dominant for determining the behavior at the macroscale. 
We note that these different categories of forces are ambiguous and not disjoint. 

Electrostatic repulsion: Due to the solid phase being e.g. negatively 
charged, the cation and anion fields at the microscopic scale in the vicinal 
fluid are neither equal nor uniform, and as a result, there is a microscopically 
varying electric field. One could solve for the electrostatic condition coupled 
with diffusion of ions (Poisson-Boltzmann equation [31, 34]), but for practical 
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situations, the complexity of the microstructure makes this a difficult task. 
The repulsion forces become significant when the increased cation/anion 
concentrations near each surface (the double layer consisting of ffist the cation 
dominant layer and then the anion dominant layer) begin to interact (double- 
layer overlap) due to the proximity of the two surfaces. These are considered to 
dominate at long-range scales. 

Van der Waals attraction: This is an attractive force acting between all 
atoms and molecules, regardless of whether they are charged or uncharged [25]. 
The current trend is to label any additional non-pressure forces not attributed to 
electrostatic forces as Van der Waals: London forces, dispersion forces, charge- 
fluctuation forces, and induced- dip ole induced-dipole forces [25]. Although some 
of these listed forces may be repulsive forces, the net Van der Waals forces are 
considered to be attractive and act on a shorter spatial scale than electorstatic 
but not as short as surface hydration forces [31]. 

DLVO theory, named after its founders, Derjaguin, Landau, Verwey, and 
Overbeek [18, 35] incorporates electrostatic double-layer forces and van der 
Waals dispersion and was developed to describe particle interactions [12]. This 
model has been criticized [:"*]. Although both forces contributing to DLVO 
(electrostatic and van der Waals) are static in nature, these two forces often 
equilibriate rapidly relative to other forces and so this assumption is appropriate 
for many systems [25]. 

3 Macroscale Forces 

Many experiments are performed at the macroscale, i.e. scale at which one 
cannot distinguish between the liquid and solid phases. Terms used at this scale 
include osmotic repulsion, surface hydration, and disjoining pressure. 

Osmotic Repulsion: Osmotic repulsion is the force that measures how 
different species interact, and is usually measured through the osmotic pressure 
experiment. For example, consider a solution (e.g. water and sugar) separated 
by a semipermeable membrane that allows water but not sugar to pass through. 
There is a difference in height and this is related to the osmotic pressure 
(technically one has to take into account the effect of the membrane but for 
the purposes of this paper we will consider it to be the difference in height). 
The osmotic pressure, vr, is the pressure that must be applied to the mixture to 
stop the influx of solvent [2, 14]. This definition holds whether one species is 
charged or not. 

We can derive an expression for the osmotic pressure. In an osmotic pressure 
experiment, the chemical potential on either side of the membrane is equal at 
equilibrium. For a component of a liquid solution which behaves as an ideal gas 
in the gaseous phase, the chemical potential is given by [2] (see also Appendix 
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where /x'^ is the mass chemical potential (energy per mass) of species j in the 
liquid phase, C'^ is the mass concentration of species j in the mixture, fip is 
the mass chemical potential of pure species j at the same temperature and 
pressure in the liquid phase, R is the universal gas constant, T is the absolute 
temperature, is the molar mass (mass of one mole of j), is the partial 
pressure of species j in the gaseous phase in equilibrium with the mixture, and 
Pm is the maximum partial pressure of species j in the gaseous phase obtained 
when in equilibrium with pure species j in the liquid phase. The activity, o}^ is 
defined as the ratio of these two partial pressures, a'^ = p^^ /pm. If the hquid 
mixture is ideal (so that Raoult's law applies) then the activity may be replaced 
with the molar concentration, 

Let's assume that on one side of the membrane the mixture is pure solvent 
(e.g. water), which we label the A^th component, and does not contain species 
j and on the other side the mixture contains species j and solvent. On the side 
of the mixture the pressure will be higher, by an amount proportional to the 
osmotic pressure, tt'^. The chemical potential of the solvent must be equal on 
both sides and we have 

RT 

^i'^^ (T, p) = /i'^ (T, p + tt'^ , C'O = /ip'" (T, p + vr'O + — ln(a'^ ) . (2) 
To evaluate ij}^{T,p + vr), we begin with the total differential 

H'' = i§rdT + J^dp. (3) 

For a pure substance, = l/p'-' where p'-' is the specific density of species j 
in the hquid phase with units of mass of j per volume of j (see Appendix A, or 
[2, 14]). Integrating at constant temperature from the state at pressure p to the 
state where pressure is p + vr we get 

/i^- (T, p + tt'O - f^l^ (T, p) = r ^ :^dR (4) 

Jp P 

Using this expression in (2b) to eliminate ij}^{T,p + tt'^) and then subtracting 
(T, p) from both sides gives: 

J-dR = -t^Ha'^). (5) 

If the density of the solvent, p''^, is constant, then we have 

1 , RT 



-TT 



piN mN 

and further if we have an ideal solution [2, 1 1], then 



ln(a'0, (6) 
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where is the molar fraction of solvent, given by moles of solvent per moles 
of mixture. If the soution is dilute, so that x'^ = 1 — x'^ where x'^ is small, then 
ln(x'^) ~ — and approximating the number of moles of N as being equal to 
the moles in the solution, we get 



where Cm is the molar concentration of j (moles of j per moles of solution), 
which is known as the van't Hoff equation. 

Surface Hydration: These are short-range bonding forces between the 
solid surface and the water that causes one to ten layers of water to be held 
tightly [25, 31]. These forces perturb the vicinal hquid, so that it behaves 
differently from its bulk-phase counterpart (water free of adsorptive forces) 
[28, 21]. Experiments by Low [ ] indicate that the macroscopic effects of 
these interactions can qualitatively completely account for many macroscopic 
experimental results. 

If one considers the solid-liquid mixture as a mixture itself, then the 
hydration forces can be thought of as the osmotic force of the solid particles. 
In fact the reverse osmotic swelling pressure experiment is presicely what was 
done by Low [ ] for montmorillonite soils. In this experiment (see Figure 
1), the liquid mixture and liquid mixture with well-layered clay minerals were 
separated by a semipermeable membrane which did not allow the clay minerals 
to penetrate, and the pressure required to keep the clay mixture from swelling 
was measured. In this case, the quantity of clay mineral was measured in terms of 
the distance separating the clay platelets (A') which represents the concentration 
of the clay mineral. In this case, the hydration pressure was determined 
experimentally to be exponentially related to the clay mineral concentration. 
If A'^ is the thickness of the clay plates then [28] 



where po is the reference (atmospheric) pressure. Equation (9) was also obtained 
via hybrid mixture theory, [l]. Note that this result is quite different from a 
pure liquid mixture in which the osmotic pressure is proportional to the log of 
the concentration, (7), but they have the same general shape - as the moisture 
goes to zero, the swelling pressure goes to infinity, and as the moisture content 
goes to 1 (A'' = 0), the sweUing pressure goes to zero. 

Disjoining pressure: is a concept traditionally used in the field of foams 
(gas-hquid dispersions) and emulsions (liquid-liquid dispersions), where the 
stability of the system relies on the stability of the thin liquid films [ ]. If 
two interfaces (in the case of foams, air-liquid and liquid-air) are separated by 
a distance h, then if h is small enough there is no portion of the interlayer (i.e. 
hquid film) which possesses the properties of the bulk fluid (see Figure 2). In 
such a case, Derjaguin and Churaev [17] state 

in mechanical equilibrium the disjoining pressure, vr (/;,), is equal to 
the difference existing between the component, Pzz of the pressure 
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Figure 1: Reverse osmotic swelling pressure experiment. 



tensor in the interlayer and pressure, Pb, set up in the bulk of the 
phase from which it has been formed by thinning out: 



7l{h) = P,, - P^ = P^ - Ps. 



(10) 



In the simplest case of a one-component liquid phase, mechanical 
equilibrium under isothermic conditions implies thermodynamic 
equilibrium. In that case the disjoining pressure is a single-valued 
function of the interlayer thickness, h,..." 

This definition has been extended so that it applies to curved surfaces by 
Kralchevsky and Ivanov [ ] . This mechanical definition is thought to be 
equivalent to the thermodynamic definition in terms of the Gibbs free energy, 
G, as [12, 19] 

7f{h) = 



dh 



T,P,A,Ni 



where the variables held fixed while taking the partial derivative include 
temperature, T, pressure, P, the area of the interface. A, and the number of 
moles of each constituent making up the thin film, iVj. 

According to Bereron [ . j] , 
the disjoining pressure is 
thought to be due to many 
forces: electrostatic double- 
layer, van der Waals disper- 
sion forces, short-range struc- 
tural forces such as hydra- 
tion, and other forces. In the 
field of thin liquid soaps, most 
treat these forces as being ad- 
ditive, although it is not clear 
that this is a valid assump- 
tion [3, 4]. More than one au- 
thor has come to the conclu- 
sion that the swelling pressure 

and average disjoining pressure are the same, e.g. [L.] (p. 282). 
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Figure 2: Cartoon illustrating variables used to 
determine disjoining pressure. 
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4 Flow in Terms of Chemical Potentials 



We begin with a formulation developed using hybrid mixture theory (HMT) 
[5, 6]. In this approach, the microscale field equations (conservation of 
mass, linear and angular momenta, energy, and electroquasistatic form of 
Maxwell's equations) are volume averaged to produce macroscopic quantities 
and equations, and then macroscopic constitutive equations are obtained by 
assuming a set of constitutive variables are a function of the same set of 
(macroscopic) independent variables and then exploiting the entropy inequality 
in the spirit of Coleman and Noll [15]. This approach has the advantage of 
developing constitutive equations directly at the macroscale, however coefficients 
in the macroscopic constitutive equations are not directly linked to microscopic 
quanitities. In what is presented here, the only geometric information retained at 
the macroscale is the volume fraction, although this approach can be expanded 
to incorporate e.g. interfacial surface density [7, 22]. In [9] the independent 
variables included 



e 



' T, v''^ S% E, Ve', VT, Vp"%... (12) 



where is the volume fraction of the liquid phase, T is temperature, a represents 
the phase (a = / for liquid and a = s for solid), p"^ is the density of the a- 
phase (mass of species j in the a phase per unit volume of the a phase), v'"'^ 
is the velocity of the liquid relative to the solid phase, £^ is the strain of the 
macroscopic (smeared out) solid phase, E is the electric field, and z°'^ is the 
fixed charge density associated with species j of phase a. 

The thermodynamic definition of liquid pressure is given by 



(13) 



where tp'' is the intensive (per unit mass) Helmholtz potential. One can either 
enforce electroneutrality with a Lagrange multiplier. A, or include an electric 
field. In the former approach, A is the streaming potential. In [6] it is shown 
that + ql.A where is the charge density and A is a streaming potential, 
is related to one third the trace of the macroscopic liquid cauchy stress tensor 
- thus the thermodynamic definition is related to what is physically measured 
[11]. 

Another pressure, the "swelling pressure", is thermodynamically defined as: 



I I I 
71 = e p 



(14) 



where is the liquid volume fraction and where the partial derivative is 
evaluated keeping the other independent variables (density, concentrations, 
temperature) fixed. It is defined so that it is a positive quantity for a swelling 
mixture. Clearly this is a macroscopic form of the thermodynamic definition of 
the disjoining pressure, and in fact, if the solid phase is structured so that it 
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does not support stress (e.g. parallel platelets), it can be shown [ ] that for a 
single component liquid, 



(15) 

G'... 



where the partial derivative is evaluated keeping the Gibbs potential (chemical 
potential) fixed, which is exactly the reverse-osmotic swelling potential 
experiment used to measure the osmotic force. Note that if the material is 
not swelling, then the energy of the hquid phase would not change with liquid 
content and the swelling pressure is zero. 

It can be shown that and tt' are related via a third thermodyanic property 
which is related to the change in Helmholtz potential with respect to volume 
keeping the mass fixed [ii]: 



I I I 9il>^ 



r (16) 



Equation (16) is mathematically exact (no assumptions), and if one converts 
to extensive variables one can show that this new quantity is the traditional 
thermodynamic definition of pressure: change in energy with respect to volume 
keeping the mass fixed. Thus the pressure in the liquid phase has two 
components: one which is the 'classical' pressure for a bulk fiuid, and the other 
the swelling pressure [11]. If the swelling pressure is zero, then the traditional 
thermodynamic pressure is the same as one third the trace of the cauchy stress 
tensor of a liquid. 

Assuming: (i) terms involving the polarization vector field are negligible, 
(ii) the gravitational term is negligible, (iii) isothermal conditions, (iv) sufficient 
moisture so that the liquid phase does not support shearing forces, (v) the charge 
associated with each species, z\ is fixed, and (v) not assuming charge neutrality, 
the resulting Darcy's law using Hybrid Mixture Theory is given by [9] 



N 



R-v^'' = _e'Vp'-7r'V£' + £'g^^;-^r'^t;'^'' (17) 

i=i 

= -£VVG" + 5^£'(pf ;^VC'^+£'g^^;-^r'^t;'-' (18) 
i=i j=i 

where ql is the charge density of the liquid phase, = ip^ — ]} j p' is the Gibbs 
potential for the liquid phase, C'^ is the mass concentration (mass of species j 
in the liquid phase per mass of liquid phase), and 1;'^'^ = v^^ — is the diffusive 
velocity. The last term involving the diffusive velocities captures the effects 
of ion hydration and relative friction between the mass-averaged velocity and 
species velocity. If the diffusive velocities (v'^'') are small then this term may 
be neglected. Note that in (17) there are no terms directly involving chemical 
potential or concentrations of species that contribute to flow. 
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Changing the concentrations of the species making up the hquid phase 
changes the pressure through (13) and through relative velocities. The 
generalized version of Pick's law for diffusion is given by [10, 31] 

vh^i = Qh . V/i"^ (19) 

where is a diffusion coefficient tensor which may be a function of volume 
fraction, temperature, and densities so that 

N 

R.v^>^ = -e^Vp^ -Ti^Ve^ + e^q[E ~^r^^Q^^ -V^i^^ (20) 
= -e'p'VG' + y e\p^f^VC'^ + e'q[E 

N 

-^r'^Q'^ ■ V/x'^-' (21) 
i=i 

From equation (20) we see that if vr' is not zero we have a pressure gradient 
threshold - i.e. gradient in the volume fraction can offset the pressure in the 
liquid pressure until vr^ Ve^ is maximum, and then further increasing the pressure 
gradient will induce flow [32, 36, 33]. An analagy between this and concentration 
gradients can be made - flow is induced by a "concentration" gradient of the 
solid phase. 

Next we express the flow equation in terms of liquid chemical potentials 
because the electro-chemical potentials are continuous between vicinal and bulk 
fluids and because it may be more useful for numerical solutions [30, 24, 23]. 

Within HMT, the chemical potential is given by [13, 10]: 

a, ^ V ) ^22) 
whereas the electrochemical potential [31] is given by 

p^". =;x"^ +^".(^ (23) 

where is the electric field potential and z"'^ is the charge density (per unit 
mass) for species j in phase a. 

We consider two cases: one in which the liquid (and bulk) phase is composed 
of only one constituent, and then a multi-constituent liquid phase. 

First consider a liquid phase which is composed of only one constituent. The 
relationship between the Gibbs potential and chemical potentials is given by 
[13]: 

N 

= ^C"^7i"^- a = l,B. (24) 
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So for a single component phase, the concentration is 1 and VC'^ = 0. Thus all 
diffusion velocities are zero, and the Gibbs potentials for the vicinal and bulk 
phases are equal up to the Lorentz term: + 5'" = where we assume the 
bulk phase fluid is charge neutral (V0^ = 0). Using the relationship between 
the Gibbs potential and the Helmholtz potential, G = ip — p/ p we get that the 
right-hand side of Darcy's equation (21), not including the hydration terms, is 
given by: 

-e^p^VG^ + e^q[E (25) 
= -eV'V^^ + a'p'/j^Vp"" - eV'^Vp^ + e'q'^E. (26) 



ipBy - - pB 

Now assume the bulk phase Helmholtz potential is only a function of density. 
Then using the thermodynamic definition of pressure, (13), the flow equation 
can be written as 

R--'' = --V|^vp- + .v^vp--.v-lvp- + .',^^;. (27) 

= -^Vp^' + e'qiE. (28) 
P 

and we see that if we write the flow equation in terms of potentials of the vicinal 
fluid as in equation (17), we have both a pressure and volume fraction potential, 
but if the flow equation is written as the bulk phase we have only a pressure 
potential. This implies that if one is measuring a vicinal pressure, there may be 
a pressure gradient threshold, but if one is measuring pressure of the bulk phase, 
there is no pressure gradient threshold. 

Now consider a multi-component liquid phase. Beginning with the right- 
hand side of Darcy equation (21) and using (24) we have: 

^ ^ ail -'V 

-^y J]V((:7'^/i'0 + $^^'(pT|^VC''^ +^'9e^-I]^'^Q'' ■ V/x'^ (29) 

j=l j=l P' j=l 

Using (22) for the thermodynamic definition of the chemical potential to 
eliminate ^t- we have (29) is equivalent to 

N N N 

-e V' C^' V/i^^ - £ V' P^' VC'j + Y ^ V' [P^' - VC'^ (30) 
j=i j=i j=i 

N 

+e^q[E - ^r'^Q^^ ■ V/i'^ (31) 
i=i 

Noting that XljLi G^^ = 1 so that the second part of the third term on the right 
side is zero and that the electrochemical potentials between the vicinal phase 
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and bulk phases are equal, we can rewrite the flow equation as: 



TV N 
TV N 

= - ^ e'p'^ V/i^^ +qfE-J2 <5'' ■ ^f^^' ■ (32) 
i=i i=i 

So in a multi-component fluid, it is the electrochemical potentials of the liquid 
phase species, Jl"^ = fi"^ + z'^^cp, that are the primary potential for fluid flow. 

To get one flnal form of this equation, consider the form of chemical potential 
for a hquid to be (see Appendix B): 

1 RT 
/i'^- = fi'^{T,po) + —ip'- Po) + — Ina^ (33) 

where pg is the speciflc density of species j in the liquid phase (mass of Ij 
per volume of Ij) and it is assumed that each component of the liquid phase 
is incompressible. Then flow equation (32) in terms of bulk variables can be 
written as: 

^ FIT ^ 

R . ^ _^pB _ y ij + qfE - V r'^ Q'^ ■ V/i^^ . (34) 

i=i i=i 

EE 
— — = 1. So in the reservoir bulk fluid 

j=i Po 

the primary driving forces are the activities (which are closely related to 
concentrations) and just as we had for a single- component fluid, the bulk phase 
pressure. Hydration of ions is an additional component whose magnitude for 
many problems has yet to be determined. 

Next we illustrate the insight that can be obtained by writing the equations 
in terms of the chemical potential. Consider Figures 3, 4, and 5, where we have 
a reverse osmotic swelling potential experimental set up where the two bulk 
phases are separated from a swelling porous material such as montmorillonite 
clay mixture by a semipermeable membrane. Across the membrane the electro- 
chemical potentials are continuous [13, 31]. Assume the chemical potential of 
the bulk phase is determined solely by the bulk phase pressure and the chemical 
potential of the liquid in the clay mixture is determined by the pressure and the 
liquid volume fraction (i.e. gradients in liquid concentrations and macroscale 
electric flelds are negligible). In Figure 3 there is no gradient in the chemical 
potential of the liquid phase, so there is no flow. The swelling pressure, vr, is 
proportional to the difference in height of the mixture and the bulk fluid. In 
Figure 4 a pressure is applied to the clay mixture. If the clay mixture is well- 
layered then (up to the hydrostatic pressure) the applied pressure is equal to 
liquid pressure, p'. However in this case the gradient in the applied pressure is 
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offset by a gradient in volume fraction, and because tlie chemical potentials in 
the bulk phase (which is in chemical equilibrium with the vicinal phase) are the 
same, there is no net flow. In 5 a pressure is applied to the left side, changing its 
chemical potential. This applied pressure causes the clay to swell more on the 
left side and because there is an overall chemical potential gradient the liquid 
flow is induced. We could recreate these same pictures for a multi-component 
liquid by keeping the pressure fixed and changing the concentrations. 



5 Comparison with Other Models 

In this section, we show how the chemical potential formulation for flow 
is a generalization of several other models, including that derived using 
homogenization (Moyne and Murad [3U]), and a mixture theoretic derivation 
of Huyghe and Janssen [23]. 



5.1 Model of Moyne and Murad 

In [30], Moyne and Murad use homogenization to upscale microcopic field 
(conservation laws and Maxwell equations) and constitutive equations to 
the macroscale. This approach provides first-order equations with precise 
expressions for coefficients in terms of solutions to the microscale equation on a 
periodic structure. The microscopic equations include: an incompressible liquid 
phase (composed of a liquid, a cation and anion) and a linear elastic solid phase; 
the conservation of momentum with the Lorentz term added {q^E where is the 
charge of the liquid (solvent) and E is the electric field); Gauss' law assuming 
polarization is negligible [ee^o'V ■ E = ql, where e is the relative dielectric 
constant of the solvent and Eq is the vacuum permittivity); the conservation 
of mass for each ion, j, dn^ /dt + V ■ j-' = with the ion flux given by 
jj — _ I (^kT)'Vji^ where Ji^ is the chemical potential of j (per molecule 
j), is the volumetric concentration of ion j, k is the Boltzmann's constant, and 
T is absolute temperature (assumed constant). In addition, electroneutrality is 
enforced. 

The solid phase is assumed to be platelet shaped (as in, e.g. montmorillonite). 
Before upscaling, a change of variables is performed to replace variables which 
may change very rapidly between the platelets (vicinal fluid) to variables that 
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are more smoothly varying. With this in mind, instead of using the chemical 
potentials of the ions in the vicinal liquid, the chemical potential of the ions 
in the bulk fluid in thermodynamic equilibrium, /i-^, is used. The apparent 
"bulk" phase fluid is a hypothetical fluid in electrochemical equilibrium with 
the vicinal fluid [j (»]. The relationship between Ji^ and /i'^ is obtained assuming 
the charged particles satisfy the Boltzman distribution, which itself assumes a 
single flat double layer. Similarly, the liquid pressure is replaced by a form of 
the "disjoining pressure": 

n = Pt-p (35) 
where p is the pressure of the vicinal fluid, and 

Pb = PlJ'w =P+ q{4>) d(p (36) 
Jo 

is the "local apparent bulk phase pressure", and where p is the density of the 
liquid, /i^ is the chemical potential of the water in the liquid phase (per unit 
mass), and is the electric potential. Note that VP;, = Vp — qeE so that Pb 
incorporates the Lorentz term. Because the relationship between p and Pb are 
assumed, it is not clear at this point whether vr is the same as the mechanical 
definition of the disjoining pressure as defined previously. After homogenizing, 
the resulting Darcy-type law is 

vl = -KV,pl-K+V,n+''-K_V,n^'' (37) 
= -pKV,./i° - K+V,n+° - K_V,n,-° (38) 

where K, K_, are second-order tensors, /i° is the chemical potential of the 
water in the liquid phase defined so that pfi^ = and a superscipt denotes the 
first term in a series expansion of orders e. In this case p^ incorporates the first- 
order approximation of the Lorentz term and the remaining terms come from 
second-order terms (fluctuations within the vicinal layer) and their relationship 
to the Boltzman distribution. 

We now show that these results are a special case of (32) Beginning with (1) 
and proceeding as we did to derive (4) we have 

/.^^(T,p^C^O = /.?(T,p) + ^lna^^- 

1 RT 
= y^J^(T,p,) + -^(p^-p,) + — lna^^ (39) 

where Ps is the standard reference pressure (1 bar) and where we assumed the 
specific densities, Pq\ are constant. At constant temperature we thus have 



V.«-=J,V.« + J^V.«>. (40) 



Let us consider that we have three species: water, j = w, cations, j = +, and 
anions, j = —■ If the solution is dilute then it is reasonable that the solvent, 
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water, will follow Raoult's law very well, so that ~ x^^ ~ 1, so that we have 
from (32) 

N 



where V is the volume of the Representative Elementary volume and n'^ is the 
number of moles of j in V. If Raoult's law applies (a^^ = x^^ , then (41) has 
the same form as (37). 

Comparing the two approaches through equations (38) and (32), we see the 
results are the same up to the definition of coefficients if we make the following 
observations/assumptions: (i) in (32) neglect the effects of hydration in the HMT 
approach, (ii) in (38) recall that the definition of Pf, incorporates the Lorentz 
term (36), (iii) in equation (32) we can assume the bulk chemical potential is 

primarily a function of the concentrations so that V/i^^ ~ l^ajVC^J , and (iv) 
the coefficients K in (38) are a function of the microscopic geometry and so are 
a function of the volume fraction. 

5.2 Model of Huyghe and Janssen 

In [23], Huyghe and Janssen use a mixture theoretic approach in a Lagrangian 
framework to develop equations subsequently used in biological applications e.g. 
[20]. 

They label their model the Quadriphasic model because they treat the 
system as consisting of four phases: cations (+), anions (-), a charged solid 
(s), and a fluid (f). Each "phase" is considered incompressible, and that 
the volume fraction of the anions and cations are negligble relative to the 
volume fraction of the solid and fluid phases. Chemical interactions are 
neglected and electroneutrality is enforced. A work energy function is assumed 
with independent variables consisting of the Green strain, Lagrangian form 
of the volume fraction of the fluid and ions, and the Lagrangian form of the 
relative velocities. The generalized Darcy law derived by exploiting the entropy 
inequality, neglecting inertial and gravitational terms is (equation (6) from [20] 
or equations (52) and (53) from [23]): 

^i^i.s = -F K [V(j9' - 7f) + n+V/I+ + n-VJi-] (42) 

where F is the deformation tensor used to convert between eulerian and 
Lagrangian frameworks, is the hydrodynamic pressure, vf is the osmotic 
pressure of the ions, /I" is the electrochemical potential incorporating the 
streaming potential, and is the volumetric concentration of phase a. Here 
the osmotic pressure, vf is deflned to be — P^, is assumed to be due to 
concentrations of cations and anions and is assumed to have a modified form of 



N 



Po 



B 



(41) 



14 



the van't Hoff equation (8), p'- — = RT(j){c^ + c~) + ttq where (j) = g\^'^Jy 
is the activity of the water or solvent, is the molar concentration of the 

water, and c' are the moles of ions per volume of fluid in porous material. 
Using n = p^ — p^ 'm. equation (42) we see that — vf is the bulk phase 

pressure. This form of the equation can be derived from (32) if one uses (33) 

for the chemical potential of the liquid phase, assumes the density of water is 

one, and neglects the hydration of ions. 

There is some question as to how well the "osmotic pressure" , p^—p^, which is 

physically the swelling pressure, can be approximated by the van't Hoff equation 

which is used for species (and not swelling) osmotic pressure. 



6 Summary 

In this paper we show that the most general way to write Darcy's law for swelling 
soils is in terms of gradients of chemical potentials, see equation (32). In this 
way one can float between using variables such as pressure and moisture content 
and electro-chemical potentials of either the vicinal fluid or a bulk fluid in 
electro-chemical equilibrium. Using this formulation tells us more easily the 
assumptions used for other models, e.g. Moyne and Murad [30] and Huyghe and 
Janssen [23]. 

One clear consequence of this formulation is determining when a pressure 
threshold gradient may exist. This is the pressure gradient that must be 
exceeded before flow is observed. A pressure gradient has been shown to exist 
when whetting a previously dry sample, [36, 33], or for swelling soils such as clay 
[32]. In this paper, we demonstrated that concentration gradients are negligible, 
if the pressure is the pressure of a reservoir in electro-chemical equilibrium with 
the swelling porous media, p^, then there is no pressure threshold gradient. 
However if the pressure is of the vicinal liquid within the porous media, p\ then 
a threshold gradient may exist - depending upon whether the swelling potential, 
TT is nonzero. See equations (20) and (28) and the discussion directly following 
them. 

This has implications in measuring the pressure - if one measures the pressure 
within a swelling porous media with a device which takes in (any) amount of 
fluid, then the fluid which is no longer affected by the presence of the charged 
solid phase becomes a bulk-phase fluid and is no longer at the same pressure as 
the fluid within the swelling porous media. Using such a device will not indicate 
a critical pressure gradient. One way that one can obtain the pressure within a 
swelling porous material is by measuring the overburden pressure - if the solid 
phase supports no portion of the stress (i.e. it is at the same pressure as the 
fluid) then the overburden pressure is the pressure in the fluid and also in the 
solid - see for example. Figures 3, 4, and 5. 
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Appendix A: The change in chemical potential with respect to 
pressure 

In this appendix we go through the calculations to show that the partial 
derivative of the chemical potential with respect to pressure while holding 
concentrations and volume fraction fixed is constant if the specific densities 
of each component is fixed. We first show this is vahd in terms of extensive 
variables and then verify the result with our definition of chemical potential in 
terms of intensive variables. 

In this section we suppress the notation for phase, a, as the definitions 
involved do not directly depend on which phase or the volume fraction of the 
phase. We assume there are N constituents making up the phase, and we define 

to be the mass fraction of component j with units (mass of j)(mass of phase). 

Extensive Variables 

Let G be the extensive Gibbs potential, G = G{T,p,M\X) where p is 
pressure, is the mass of species j, and X is any other variable upon which 
the Gibbs potential depends, such as the volume of the porous media. We note 
that normally we write G as a function of the number of moles of species j, N\ 
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but the ration of and A^-' is the molecular weight of j (with units of mass 
of j per mole of j), and since the molecular weight is a constant this does not 
change the following results. 

The thermodynamic definition of chemical potential in units of energy per 
unit mass is given by 



dG 



(A.l) 



T,p,X 



Before deriving a Maxwell relation we use the total differential to determine 



dG = d{U-TS + pV) 

= dU -TdS - SdT + pdV + Vdp 



N 



31 J f)U F)JJ r)U 

dS + —dV + V T^dM^ + —dX - TdS - SdT + pdV + Vdp 



dS 



dX 



N 



BU dU 

TdS -pdV + J2 ^TTjdM^ + —dX - TdS - SdT + pdV + Vdp 



dX 



N 



dU BU 



dX 



(A.2) 



where we used T = and |^ = —p [13]. 

Now taking the partial of both sides with respect to p keeping the appropriate 
variables fixed we have: 



dG_ 

dp 



V. 



(A.3) 



Now let's assume that the function G is smooth enough so that mixed partials 
commute. Using (A.3) we have 



d^G 



d'^G 



dp dM^ dM^ dp 
d f dG\ d fdG\ 



dp \dM^ J 



- 



dp 



dM^ \dp J 
dV 



T,M'',X 



dM^ 



(A.4) 



T,p,M>={k^j),X 



Equation (A.4) tells us that the partial of the chemical potential with respect 
to pressure is related to how the volume changes with the quantity of j. 

Recall that our goal is to show that if the intensive densities of the species 
(so the mass of species j per volume of species j) are constant then so is 
The units of fP are mass of j per volume of mixture. So let's define the intensive 
density to be pg; which has units of mass of species j per volume of species j. 
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We use a subscipt to emphasize the fact that it is not equal to fP and that in 
what follows we consider to be a constant. 

Let V = + V'^ + • ■ ■ + V'^ be the volume of the mixture, where is the 
volume of species j. Then we have = PqVK With this, (A. 4) gives: 



dp 



dV 



dMi 

N 



N 

E 

k=l 



dV'' 



dMi 



T,p,M'{¥j),X 



dV^ 



N 



l_ dV^ 
2-^ J 



k=i Po 

1 



dV^ 



T,p,Vi{l^j),X 



T,py\i^j),x 



P'o' 



(A.5) 



So if the specific densities for every component is fixed then the dependence of 
the chemical potential upon pressure is linear. Note that if only one species is 
incompressible, then we would need the additional assumption that the density 
of species j is independent of the quantity (volume) af all species (including j), 
in order for (A.5) to hold. 

Intensive Variables 

We now go through the same argument in terms of intensive variables. Since 
the definition of chemical potential is relatively new, [10], we go through the 
calculations in detail to verify the same result holds. The chemical potential as 
defined in this paper in terms of the Helmholtz potential, (22), is written with 
assumed independent variables, T, pP , X, where X could be any other variable. 
We would like to determine the definition of chemical potential in terms of the 
Gibbs potential, g{T, p, C^). To do this we first look at the definition of chemical 
potential in terms oi ifj = '4'{T, p, , X). 

Claim: 



N 



i=i 



T,p\l^j),X 



2 ^ 

dp 



(A.6) 



T,C^{l=l,..,N-l),X 



To show this result we begin with the equivalencies of the total differential 
of the Helmholtz potentials: 

v^(T,p,c^x) = v>(T,p^x) 

#(T, p, X) = #(T, p^ X) 



N 



MKT, p,C',X) = ^dT + j: ^dp" + ^dX. 

k=l p 
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Now taking the partial with respect to p on both sides keeping the 
concentrations (and T and X) fixed we have: 



dip 
dp 



N 



dip dp^ 



T,CHj=l,...,N-l),X 



k 

N-1 



^ dp'' dp 



^ dip diC'p) 



k 



^ dp'' dp 
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iV 



d 



dip 



k=l 
N 



dp'' 



dp 



N 



^ dp''' 

k=l ^ 



:^-Y.tiC')p 



dp 



Multiplying both sides by p^ and using the fact that C''p = p'' we get (A. 6). 



Claim: The chemical potential in terms of the Helmholtz potential, 
is given by 



p^ 



d{pip) 



dp^ 



N-l 
P ^ ^k 



^ k=i 



dip 



dC' 



+ 



dip 



T,p,C'{l^k),X 



dC^ 



j = l,...,N (A.7) 



where S^'^ is one if j = and zero otherwise. 

We begin as we did in the^ previous claim by equating the two functions 
of Helmholtz potential, ip = ip and looking at the total differential. We will 
then use the thermodynamic definition of chemical potential given by (22), 

dpi 



P' 

dip 



p'={kytj) 



dip 

dip 
df 



p,C'=,X 



dp 



T,C^,X 



k=l 
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dC^ 



T,p,C'{l^k),X 



+ dX 



dX. 



Taking the partial derivative of both sides with respect to fP we have 



dip 
dp^ 



dip dp 
dp dp^ 



N-l 



dip dC' 



^ dC' dpi 



T,pi{l^j),X 

We now need to evaluate the terms dp/{dp') and dC'' / {dp^): 

N 



T,p,C'' 



(A.8) 
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Also: 



j = k: 
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(A.IO) 



where we used (A. 6) and is the result of this claim. In the above we note that 
ip = ip a.s these represent the same quantities and we can choose the functional 
form of the Helmholtz potential. 
Also note that we have 



7 N 
P^ - P 



dip 



dC^ 



(A.ll) 



T,p,C^{lytk),X 



which is the relationship derived using an exploitation of the entropy inequality 
in [10]. 

Claim: The chemical potential in terms of the Gibbs potential is given 
by 



N-l 



P- 9 



dg 



k=l 



dC^ 



T,p,C'{ljtk),X 



99 
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6^ 



j = l,...,N (A.12) 



where 5i'^ is one if j = and zero otherwise. 

We first derive the intensive equivalent to (A. 3) by beginning with the 
relationship between the Gibbs potential and the Helmholtz potential, 

9 = '^ + -, 
P 
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where g = g[T,p, , X) for /c = 1, . . . , — 1. Taking the total differential of 
both sides: 



dg 
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dp 
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(A.13) 



where we used (A. 6) to cancel two terms in the last step. 

To get the equivalent of (A. 3) take the partial with respect to p on both 
sides and we have 



dg_ 
dp 



T,Ci(j=l,...,N-l),X 
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(A.14) 



which is consistent with (A. 3) in the sense that if we divide both sides of (A. 3) 
by the total mass (and the total mass is fixed) we get (A.14). This remark just 
shows consistency. 

To get the chemical potential in terms of the Gibbs potential begin with 
(A.13), 



dg 



dC^ 
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T,p,C'^{kjtj),X 
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(A.15) 



T,p,CHkjtj),X 



and use the result from the previous claim, (A. 7), to get (A. 12). 

Two checks can be made on this result. If there is only one component 
(A^ = 1) then the chemical potential of the phase should be the Gibbs potential, 
and the sum of the weighted chemical potentials should be the Gibbs potential: 

N 

Y^0p^=g. 

i=i 

A few algebraic steps shows that both of these results hold. 
Further, letting j = in (A. 12), we get 
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Claim: If each component of the phase is incompressible, then the 
partial derivative of the chemical potentail with respect to pressure 
is constant. 

We adapt the notation from the extensive results and let Po be the intrinsic 
density of component j (mass of j with respect to volume of j). For this claim 
we assume that Po is constant for j = 1, . . . , A^. 
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We begin by showing a preliminary results using = PqV^ where is the 
volume of j per unit mass of the phase material, and is not 1/ which has units 
of mass of j per unit volume of phase material. Let v = 1/ p = XljLi '^^ ■ Then 
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Now let's determine the partial derivative with respect to /i^ first. Beginning 
with (A. 16) and using (A. 14) we have 
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where we used result (A.17). 

Now the rest is easy if we begin with (A.ll) and (A. 15): fi^ 



(A.18) 
fx^ + dg/{dC^) 
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for J = 1,...,A^-1: 
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(A.19) 



where we used (A. 14) and (A. 18) in going from line 2 to line 3. and (A. 17) to 
go from line 3 to line 4. 

And so we see that if the specific densities are constant (i.e. do not change 
too much with the temperature and concentration fluctations of the particular 
problem being considered), then the chemical potential changes linearly with 
the total pressure. This result is generally used for liquids and not for gasses. 

Appendix B: Background Material on Chemical Potential 

This section contains material found in a standard textbook on physical 
chemistry [2, 14]. It is presented here for easy reference. 

The chemical potential has three defining properties (1) it is a scalar quantity 
representing the energy change as the quantity of species is changed (partial 
derivative of a potential with respect to quantity), (2) is a quantity which is 
equal in two different phases at equilibrium, and (3) is the generalized driving 
force for diffusion. 

For a pure substance in a single phase, the chemical potential is equal to the 
Gibbs potential (per unit mole), G. We first determine the Gibbs potential for 
a single compenent, ideal gas that satisifies pV = nRT where V is the volume, 
p is the pressure, n is the number of moles, R is the gas constant, and T is the 
absolute temperature. In this case, the Gibbs potential is only a function of 
temperature and pressure: G = G(T,p) and [13, 2] 



To determine how the Gibbs potential depends upon pressure, integrate the 
above relationship from a reference state Go(Tq,pq) (where Pq is the standard 
pressure) to a second state at constant temperature and number of moles, 
Go(Tq,p), and using the ideal gas relationship p = RT/vm (where Vm is the 



dG 




-SdT + VdP. 
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molar specific volume with units of volume per mole) we have 



G{T,p) = Go{T,po)+ fvmdP 

— r RT 

= Go(r,po)+ / —dP 

= Gq{T,po) + RT\yi 



PO 

' p 



Po 

If we do not have an ideal gas then we replace the pressure by an effective 
pressure, called the fugacity, f, and we have 



G{T,p) = Go{T,po)+ fvmdP 

Jpo 



= G^{T,p^) + RT\n(^\, 

\PoJ 

and in fact, this is the definition of fugacity. 

For a mixture of gases, define the partial pressure of species j to be — x^p 
where is the molar fraction of species j (moles of j per moles of mixture) . For 
a mixture of perfect gases (each gas ideal and no interactions between species) , 
the partial pressure of species j would actually be the pressure of species j if 
no other species were present (Dalton's law). Using the definition of partial 
pressure we have, for a component of a perfect mixture of gases: 



77* (T, p, x^) = Jll^ (T, p,) +RT In (^^y 



(A.l) 



where JI^^ is the chemical potential of species j in units of energy per mole, and 
Po is the standard pressure (which is 1 if pressure is measured in bars). Using 
Dalton's law, we have 

ll'^{T,p,x^) = il',^{T,po) + RT\n(^ 

\ Po 

= ]li{T,po) + RTln (^^^ + RThixK (A.2) 

Since Xj is always between and 1, we have that the last term is always negative 
and so the chemical potential of a component in a mixture is always less than 
the chemical potential of a pure substance. We should note that the above 
result is only true for a mixture of ideal gases. This model breaks down if for 
example, there are chemical reactions, the pressure is high, or there are strong 
intermolecular forces between (or among) the different species. Water vapor 
in the atmosphere is usually treated as an ideal gas, with error in e.g. density 
calculations of less than 0.2% (http:/ /en. wikipedia.org/wiki/Density_of_air). At 
high pressures the perfect mixture assumption breaks down. 
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For a mixture of nonideal gases, the partial pressure must be replaced by the 
fugacity and we have 

71* (T, p, x^) = Jll' {T, Po) + RT In 

\Po 

For a liquid, the chemical potential is determined by using the fact that the 
chemical potentials of one species in two phases are equal at equilibrium. Thus 
the chemical potential of a pure liquid of an ideal component, j, (the component 
behaves as an ideal gas in the gaseous state) is: 

Jl'^ = -fll^ = Jll^ (T, Po) + RT In (^^^ , (A.3) 

where /I* is the chemical potential of species j in the gaseous state in equilibrium 
with the liquid state, JIo is the chemical potential of species j in the gaseous 
state at the same temperature but at standard pressure (1 bar), is the partial 
pressure of j in equilibrium with pure j in the hquid phase, and po is the standard 
pressure. Now suppose we have a liquid mixture of ideal species. Let be the 
partial pressure of species j in the gas phase. Then the chemical potential is: 

71'^ = 71^^ (T, ) = 71^^ (T, po) + RT In (^] = 71'^ (T, p) + RT In 



= Jl'^iT,p) + RT\n{a^), (A.4) 

where we eliminated JI^^ using (A.3) and = p^^ /pp is the activity of component 
j. We note that for water, the activity is the relative humidity divided by 100 
(i.e. relative humidity not in percent form). Since the partial pressure for a 
species in a mixture is usually less than the partial pressure in the pure state 
(i.e. in the case of water, the relative humidity is between and 1), we see 
that the chemical potential for a species in a mixture is generally lower than 
the chemical potential of a pure species. The above equation holds whether the 
liquid solution is ideal or not. 

Now if we have an ideal liquid solution, then Raoult's law is satisfied (this is 
the definition of an ideal solution): 

where p^^ is the partial pressure of species j in the gas phase in equilibrium 
with the ideal solution, is the molar volume fraction in the liquid, and pp^ of 
the partial pressure of species j in equilibrium with pure liquid j. Using this 
relationship we have 

Jt"^ =Jl'o'{T,po) + RT\nx^, 

where again, x^ is the molar fraction of species j in the ideal solution. An ideal 
solution is one in which the liquid molecules in the mixture interact with all 
other species the same, i.e. there is no difference between how species i interacts 
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with i and how species i interacts with j. This is a much stronger assumption 
than assuming that the species behaves as an ideal gas in the gaseous phase. 
Raoult's law is known to hold for a solvent when it is nearly pure. For real 
solutions where species j is at low concentration, instead of p*' = x^pp' we 
replace pp^ with an empirical constant (measured) so that p^^ = x^K^ and this 
is referred to as Henry's law. 

Returning to (A. 4) we want to use our knowledge from Appendix A to express 
the chemical potential of a liquid in terms of pressure. Using the total differential 
and that /l'-' = m-'/i'^ we have: 

Pi ^ 

where no assumptions have been made. Integrating both sides from standard 
pressure to pressure, p, at a constant temperature we have 

Ji'^{T,p)--p^^{T,po) = — (p-po) (A.5) 

Po 

if p'^ is constant over the pressure range Pq to p. Combining this with (A. 4) we 
have 

Jl'^iT,p,x^)=-fl'^iT,po) + —ip-po) + RTln{a^). (A.6) 

Po 

So let us compare the difference between chemical potentials in the liquid 
and gas phases, (A. 2) and (A.6): 

(T, X* ) = JlUt, Po) + RT\n(^] + RT In . (A.7) 

\PoJ 

Ti';{T,p\x'^) = Jl'^{T,po) + —{p'-po) + RT\naK (A.8) 

Po 

Let's consider the case where j is water, the gas phase is atmospheric air, and 
water is the primary component of the liquid phase. In this case the ratio 
of pressures in the gas phase is close to 1 (atmospheric pressure is close to 1 
bar) and so the pressure term drops. For the liquid phase in which water is 
the primary component, the activity is approximately the molar concentration 
which is 1, and the term involving the activity is negligible. Thus we have 

p^{T,p3,x^^) ^ Jli{T,po) + RTlnx^'. (A.9) 

l2'^{T,p\x'^) ^ Ti'^{T,po) + —{p' -Po) (A.IO) 

Po' 

that is, the chemical potential of water in the gas phase is primarily determined 
by the concentration, and the chemical potential of water in the liquid phase is 
primarily determined by pressure. 



28 



If however there is something other than water in the hquid phase, then 
the chemical potential of a component in the liquid phase is determined by the 
pressure and the relative humidity in equilibrium with the liquid phase. 



Appendix C: Nomenclature 

In general, a superscript Greek letter indicates a macroscale quantity from 
that phase. Superscript minuscules indicate the constituent, so that, e.g. v^, is 
the macroscopic velocity of constituent j in the a-phase. Subscript b refers to 
the quantity in the bulk or reservoir phase in electrochemical equilibrium with 
the vicinal fluid. 

a'-'' activity of species j defined to be the ratio of pressures p^^ /pp, [-] (2) 
A area [length^] (11) 

mass concentration, p^'^/p", [-] (1) 
E electric field [force/charge] (12) 
£^ macroscale strain tensor of solid phase, [-] (12) 
F deformation tensor [-] (42) 

G° Gibbs free energy, [energy of a-phase/mass a] (11), (24) 
j-' ion flux [length/time] Section 5.1 
molar mass [mass / mole of j] (1) 

volumetric ion concentration [volume ion j'/volume of solvent] Section 5.1 
iVj number of moles of species i [moles] (11) 
p pressure [force/area] (1) 

classical pressure (1/3 trace of cauchy stress tensor at rest), [force/length^] 

partial pressure of species j in the gas phase (1) 
Pp maximum partial pressure of species j in the gas phase (1) 

charge density of phase a, ^,=iP"^-2"^ [charge a/volume a] (17) 

coefficient for generalized Pick's law (19) 
r'^' coefficient for capturing ion hydration effects (17) 
R Universal gas constant, [Force-length/(degree-moles)] (1) 
t time 

T absolute temperature [degree] (1) 

v'^i velocity of jth constituent in phase a, (18) [length/time] 
■y" mass-averaged velocity of phase a, '^j^^C°'^v°'^ , [length/time] (12) 
diffusive velocity, i;"^ — [length/time] (12) 

velocity of liquid relative to solid phase, — v'^, [length/time] (12) 
x'-'' molar fraction of j*^ species in phase /, [-] (7) 

z"^' fixed charge density associated with aj, [ charge a^/mass aj] (12), (23) 
£° volume fraction of a-phase in Representative Elementary 

Volume (REV), \6Vc,\/\6V\, [-] (12) 
£ dielectric constant of solvent [-] Section 5.1 
£o vacuum permittivity [charge^/Force-Length^] Section 5.1 
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A* thickness of montmorillonite clay mineral, [length] (9) 

A' thickness of vicinal liquid of montmorrilonite clay, [length] (9) 

/i"^ chemical potential of jth constituent in phase a 

[energy aj mass a^] (1), (22) 
/i"^ electro-chemical potential of jth constituent in phase a 

[energy aj mass Oj] (23) 
/Z"^ chemical potential of jth constituent in phase a 

[energy aj mole in Appendix B (23) 

[energy aj molecule in Section 5.1 
v}^ osmotic pressure, [force/area] (2) 
7r° swelling pressure, [force/area] (14) 

electric field potential, E = —'Vcf), [force/charge- length] (23) 
Helmholtz free energy density of a-phase, 
e° — T?7" [energy a / mass a] (13) 

density of j^^ constituent in phase a, C"^ p", [mass aj/ volume a] (12) 
p"^ specific mass density of phase a, [mass aj/ volume aj] (4) 

averaged mass density of phase a, ^j=iP"-' , [mass a/ volume a] 
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